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Abstract 

We study the contribution of instantons to the binding and the mass of the lightest scalar glueball, 
in gluondynamics. We show that the short-range correlations introduced by such non-perturbative 
vacuum fluctuations are sufficient to give raise to a scalar glueball bound-state, with mass in good 
agreement with the results of recent lattice calculations. 
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I. INTRODUCTION 



The question whether the hadron spectrum contains glueballs and hybrids has received 
much attention, since the establishment of QCD. From the experimental point of view, such 
exotic resonances are believed to be most effectively excited in photo-production experi- 
ments. However, the available data do not allow to unambiguously identify the gluonic 
exitations, among the ordinary \qq) states [l|. In order to help clarify the ambiguous iden- 
tifications below cc threshold, dedicated spectroscopic experiments with linearly polarized 
photon beams will be performed by the GlueX Collaboration at Jefferson Laboratory 
The PANDA experiment at GSI [3| has been designed to discover glueballs up to masses of 
4.5 GeV and will be fully operational by 2012. 

From the theoretical point of view, the problem of understanding the gluon content of 
hadronic wave-functions is deeply connected with the problem of identifying the dominant 
non-perturbative vacuum gauge configurations, responsible for the dynamical quark-gluon 
and gluon-gluon correlations, at intermediate and large distances. In particular, the spec- 
trum of gluondynamics encodes the information about the vacuum gauge configurations 
responsible for the non perturbative gluon-gluon interaction. 

Gluondynamics has been extensively investigated with quenched Lattice QCD simula- 
tions 0, [5]. These studies have established that the lightest glueball is a J PC = ++ state, 
with mass in the range 1.500 — 1.750 GeV. In particular, Chen et a/.[4J recently reportec 
mass of M ++ = 1.710 GeV, obtained with anisotropic lattices, while Meyer and Teper 
used much finer isotropic lattices, and reported a value M ++ = 1.475 GeV. Both such cal- 
culations are affected by a ~ 2% statistical error and a ~ 3% systematic error, associated 
with the determination of the physical lattice spacing. 

It has also been established that the physical properties of this state are remarkably dif- 
ferent from those of all other glueballs. First of all, it is significantly lighter, since tensor and 
pseudo-scalar glueballs have masses M 2 ++ — 1.4 M ++ and M -+ ~ 1.5 — 1.8 M ++, respec- 
tively. In addition, the scalar glueball is an unusually compact object, of size 1 ~ 0.2 fm jf|, 
while the 2 ++ and h states have much larger sizes > 0.8 fm, i.e. comparable with that of 
the typical hadrons. It is interesting to note that in the MIT Bag Model, in which all states 



a 
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1 The size of the glueball can be estimated, e.g. by measuring the exponential decay of Bethe-Salpeter 
wave-function in Euclidean space 
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are bound exclusively by confining forces, the scalar and tensor glueballs are predicted to 
have comparable masses of ~ 1 GeV and comparable sizes of ~ 0.8 fm. 

Given the very small size of the scalar glueball, one expects that the dominant forces re- 
sponsible for its binding must be provided by short-sized vacuum fluctuations, of correlation 
length ~ 0.2 fm. Consequently, typical energy-momentum transfers inside a scalar glueball 
occur at the chiral symmetry breaking scale A x ~ 1 GeV, rather than at the confinement 
scale Aqcd ~ 250MeV. 

It has long been argued that instantons of size ~ 0.3 fm represent the main vacuum gauge 
field configurations responsible for the non-perturbative dynamics at the chiral scale [7, 



Several recent lattice studies [9|, llOj have provided strong evidence in support of such an 
hypothesis. Hence, it is natural to raise the question whether the correlations generated by 
such semi-classical vacuum fluctuations are sufficiently strong to bind the lightest glueball 
state. 

The role of instanton-induced interactions in light hadrons has been investigated in a 
number of papers. It was shown that they giv e raise to a nucleon and a pion with realistic 



masses 
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12j and electroweak structure 



161 ]. In addition, the dependence of the 



mass of these hadrons on the quark mass agrees extremely well with the predictions of chiral 
perturbation theory and with the available results of lattice QCD simulations [14]. It was 
also shown that lowest-lying vector and axial-vector meson resonances can be bound by 
instanton-induced chiral forces. In these resonances, however, instanton effects are much 
weaker than in the nucleon and pion, as they are suppressed by a power of the instanton 
ensemble's diluteness, k ~ 0.1. As a result, the mass of the resonances turns out to be about 
30% heavier than the experimental value. This fact can be interpreted as a signature of the 
fact that other sources of non-perturbative correlations — presumably confining forces — 
begin to play an important role in these systems. We argue that a similar situation is 
realized in the glueball sector: the lightest scalar is bound predominantly by instantons, 
while the structure of the spectrum of its resonances is presumably shaped by other gauge 
configurations responsible for color confinement. 



The instanton contribution to the lightest 
work of operator product expansion (OPE) 



due 
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Is has been investigated in the frame- 



20| or by computing directly point- 



to-point correlators at short and intermediated distances, in the Random Instanton Liquid 
Model and in the Interacting Instanton Liquid Model [21]. Such studies allowed to establish 
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that, at least at the qualitative level, the ordering of the states in the light glueball spectrum 
fits well in the instanton picture: in fact, the O(k) instanton interaction is attractive in the 
scalar channel, is suppressed in the tensor channel and is repulsive in the pseudo-scalar 
channel. Remarkably, the same dynamical mechanism can explain the ordering of the light- 
est qq excitations. In fact, the 0(k) instanton correlations are attractive in the pion, are 
suppressed in the p-meson and are repulsive in the ?/-meson. Finally, the same argument 
applies also to the lowest-lying sector of the light baryon spectrum, where leading instanton 
effects are attractive in the nucleon, but suppressed in the A-isobar. It should be stressed 
that the short-range interactions induced by instantons in these systems can at most provide 
binding, but do not confine. 

At a quantitative level, the most sophisticated instanton-improved sum-rule calculation 
available [^(J leads to a scalar mass ms = 1.25 ± 0.20 GeV, while numerical simulations of 
point-to-point correlators in Interacting Instanton Liquid Model give ms — 1.5 GeV 21]. 
The fairly good agreement between these calculations may be interpreted as evidence in favor 
of the existence of scalar glueball in the instanton vacuum. On the other hand, it should 
be stressed that both calculations rely on the assumption that the spectral function in the 
scalar channel contains a glueball bound-state, separated by a "perturbative" continuum. 
The existence of a pole in the Green's function is postulated, and its location is determined 
from a fit of the calculated correlators, in either Borel or coordinate space. 

In this work, we complement these studies using a different approach. We address the 
question of the existence of a scalar glueball without assuming that the corresponding two- 
point correlation function contains a pole. To do so, we perform an analysis which is quite 
similar to the one used to extract the mass in lattice simulations. However, rather than 
looking at the long Euclidean time behavior of the zero-momentum projected correlation 
functions, we focus on finite momentum states and we perform a calculation in a finite 
range of Euclidean times. This choice allows to reduce the model dependence associated to 
the unknown structure of the instanton-intanton correlations. This way, we find evidence 
for a scalar glueball bound state, whose mass is in good agreement with the results of recent 
lattice calculations. 

The paper is organized as follows. In section [III we describe the method used to infer 
the existence of a glueball and to determine its mass. In sections III 1 1 and IIVI we present 
our calculation. Section |V] contains some general considerations on the structure of gluonic 
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correlations responsible for the binding of the glueball, in this model. Finally, the main 
results and conclusions are summarized in IVII 



II. SCALAR TWO-POINT CORRELATION FUNCTION AND THE 
MOMENTUM-DEPENDENT EFFECTIVE MASS 

In quantum field theory, the information about the spectrum is encoded in the two-point 
correlation functions. In our specific case, we are interested in the scalar correlation function 



where x = (r, x) is a space-time point in the Euclidean space and S(x) is proportional to 
the action density operator 



which excites states with J PC = ++ quantum numbers. 

Using the spectral representation, it is immediate to show that, in the limit of large Eu- 
clidean separation \x\ — > oo, the point-to-point correlator (CQ) converges to the square of the 
gluon condensate. At large-but-finite space-time distances, Hs(x) receives contribution also 
from the lowest-lying vacuum excitation, with ++ quantum numbers. Lattice simulations 
have shown that, in gluon-dynamics, the lightest vacuum excitation is a glueball. In QCD, 
if a scalar glueball exists it must be unstable and decay into a two-pion state. In general, 
in order to extract information about the scalar glueball from the large-distance behavior of 
the correlation function (CQ) one has to subtract the gluon condensate contribution and as- 
sume that the purely gluonic operator S(x) couples weakly to hadron excitations containing 
quarks. 

For the purpose of identifying a signature of the existence of a glueball bound-state in 
the spectrum, it is convenient to introduce the notion of effective energy. We consider the 
momentum-projected correlation function, 



The effective energy E e ff(r, p) is defined as the logarithmic derivative of Gs(r, p), 




(1) 




(2) 




(3) 




(4) 
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If the spectrum contains a scalar glueball bound-state of mass M ++ , then at large Eu- 
clidean times the effective energy E e ff(r,p) must become independent on r and converge 
to the kinetic energy of such a state, propagating with momentum p: 



lim%(T,p) = Jp 2 + M 2 ++ . (5) 

T— >0O V 

In lattice QCD simulations it is convenient to project onto states with zero momentum, 
p = 0, since projections on finite momentum states are usually more noisy. This way, one 
can directly read-off the mass of the lowest state from the plateau in the effective-energy, at 
large Euclidean times. 

However, one price to pay for such a choice is that one needs to subtract the contribution 
of the gluon condensate from the scalar two-point correlator. If the fluctuations of the 
correlator are very large, or if the simulation box is not sufficiently large, the result of such 
a subtraction may very noisy or affected by systematic errors. In this work we choose to 
adopt a more general standpoint and consider the projection of the point-to-point Green's 
function to some finite momentum p. Hence, we define a momentum- dependent effective 
mass, 



M eff (T,p) = ^E* eff (T,p)-p2. (6) 

If the lowest scalar excitation in the spectrum is a single-particle bound-state, then in the 
large Euclidean time limit M e jj(r, p) must stop depending on r and on p and converge to 
the glueball's mass: 

lim M e// (r,p) =M 0++ . (7) 

T^OO 

We conclude this section by discussing how the range of Euclidean times where the lowest 
excited state contribution becomes dominant depends on the projection momentum p. The 
spectral reppresentation of the momentum-projected correlator ([3]) leads to the expression 

G(T,p) = 5>„(p)e-^, 

n 

where E n (p) is the energy of the n-th eigenstate of momentum p and c n (p) is the corre- 
sponding overlap with the scalar operator. 

The momentum dependence of the couplings c n (p) is completely determined by Lorentz 
symmetry. For example, the coupling of the scalar operator to the scalar glueball state reads 

(0\S(0)\G o++ (p)) = . f° + " 2 (0|g(0)|G o++ (0)>. (8) 
yp 2 + M 2 ++ 



Let us focus on large Euclidean times for which only the two lowest vacuum excitations 
contribute. If these are single particle states, in the zero-momentum projection we have 
E% = mi an d E 2 = m 2 . Hence, the time required to filter the lowest state scales with the 
inverse mass gap r ~ — - — . 

At large momentum the filtering of the lowest state becomes more difficult, since E% = 
\/ m\ + p 2 fa E 2 = \Jm 2 l + p 2 . Consequently, the contributions of the excitations require 
longer times to die out. 

We now estimate the dependence of such a time scale on the momentum p. Let Gi(t, p) 
and G 2 (r,p) be the contribution of the lowest two excitations to the momentum projected 
correlator, i.e. 2 

G(t, p) = Cl (p)e-^(^ + c 2 (p)e-^ =: G 1 (t, p) + G 2 (r, p). 

The lowest excitation dominates if 

G 1 (r,p)>f-G 2 (r,p), 

where / is some large parameter, / 3> 1. This will be the case for times larger than some 
given time r , defined as 

logCfg) 

r ° E 2 ( P )-E 1 ( P y 

For ultra-relativistic momenta p ^> mi,m 2 , this leads to tq oc p. In Fig. [T], we show the 
expected behaviour of the corresponding effective mass at different momenta, when several 
boundstates are present in the spectrum. One can clearly see that at higher momenta, the 
plateau which can be associated to the lowest state is shifted to larger Euclidean times. 

III. SCALAR TWO-POINT CORRELATOR IN THE INSTANTON MODEL 

The classical field of an individual instanton or antiinstanton is specified by a set of 
collective coordinate Q = (z, p, U) where z is the position of the pseudo-particle, p is its size 
and U is a SU(3) matrix which selects the orientation in color space. 

Let A^(x; Q\, Qn=n + +N-) be a back-ground field, constructed from the classical fields 
of N + instantons and iV_ anti-instantons, according to some prescription. For example, the 

2 Note that we are neglecting the non-exponential dependence of the coefficients c\ and C2 on the momentum 
P 
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Effective mass plot at different momenta 

> 
o 
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FIG. 1: Structure of the momentum-dependent effective mass plot at different momenta, for a 
spectrum consisting three boundstates of mass 1 GeV, 2 GeV and 4 GeV. 

simplest choice is the so-called sum ansatz, in which amounts to adding up the contribution 
of all individual fields: 

N 

Af(x) = Y, 4^°*)> ( 9 ) 

where A 1 (x; f^) is the classical field of the z-th instanton. 

The Instanton Model for gluon-dynamics is defined by replacing the path-integral over the 
gauge field configurations by a sum over the configurations of a grand-canonical statistical 
ensemble of singular gauge instantons and antiinstantons: 

N 

Z ILM = JjIaFI / 11^ e-*~[<< ni --°*>] e-^>-^) (10) 

N+,N- +■ -■ J i 

Sym[A^'(£Ii, •••fijv)] is the classical Yang-Mills action, evaluated in the classical background 
field A c J-(p,x, O/v). F(Qi, . . . is a function which formally represents the result of 
performing the functional integral over the fluctuations around the background fields, in the 
original path integral. Based on the analogy with the classical statistical mechanics, this 
term can be interpreted as the "potential of mean-force" between the pseudo-particles, due 
to all the degrees of freedom in the original system, which have been integrated out (i.e. the 
quantum gauge field configurations). Clearly, if one could calculate F(fti, Q^) from first 
principles, then Eq.f llOp would represent an exact re-parametrization of the path- integral of 
gluon-dynamics, which would not stand on the semi-classical approximation. 

In general, both the classical action SymIA^} and the quantum effective action 
F(Qi, Q N ) introduce many-body correlations between the pseudo-particles. If the en- 
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p=0 GeV 
p=2 GeV 
p=4 GeV 



semble is sufficiently dilute, it is possible to neglect three- and higher- body correlations and 
re-write the sum as: 

Sym^, ...,Q n )+ F(Q u ...,n N )~Y, G i(^) + Yl G ^ ( n ) 

i i<j 

where the functions G\ and G 2 completely specify the model. In the appendix, we illustrate 
the choices of G\ and G2 which lead to the Interacting Instanton Liquid Model. 

Let us now discuss the evaluation of the two-point scalar correlator ([3]), in the Instanton 
Liquid Model. The first step consists in writing the gluon field operator in the operators as 
a sum of a classical instanton background and a fluctuation field, 

A,(x) = A™ t ix) + Al uct ix). (12) 

At distances much smaller than the typical instanton size, the scalar correlator is dominated 
by the contribution of the fluctuation field, since the instanton contribution is suppressed by 
the finite size of the pseudo-particles. On the other hand, at large distances, the correlator 
is expected to be dominated by the background fields, which are non-perturbative. We 
then approximate the complete scalar correlator as a sum of two terms: a short- distance 
part, estimated from free gluon propagation, and a long-distance, generated by the classical 
instanton fields [221 ]: 

n 5 (x) = n| ee (x) + n^ st '(x), (13) 

where the short-distance part is readily computed and gives 

and g is the running coupling constant. 

-Qmst.^-j - g correlator in the classical multi-instanton background field, is evaluated 
in the canonical ensemble of N + instantons and iV_ anti-instantons and reads 

1 1 1 r N 

wtix) = — ^-y— J jldnje-^Gim-^G^M c^d) (15) 

where S(x, A^) is the scalar operator (j2J) evaluated in the classical back-ground field, and 
Zj LM is the canonical partition function, 

Z! LM = [ fld^e-**™-*™^™ (16) 

+• • J j=l 



We now undertake one more approximation and assume that the scalar operator ([2]) can 
be written as the sum of the contribution coming from the individual pseudoparticles: 

S(x;Cl 1 ,...,Cl N ) =} j s(x; O fc ), (17) 
k 

where s(x; flk) is the action density evaluated at x, in the field of the k — th speudo-particle 
and reads, 



s(x; tt k = (z k , p k , U k )) = xo fc on4 - (18) 



Note that this approximation is motivated by the fact that the action density decreases very 
rapidly away from the instanton — according toal/i 8 power-law — . Hence, the ansatz ffT7|) 
becomes accurate if the system is dilute. 

With such a choice, the instanton contribution (Tl5|) can be decomposed in two parts 

1 i i r N 



3=1 



X 

Lm 



^2s(x-,Qi) s(0;fl m ) =U 1 (x) + U 2 (x) (19) 



lli(x) represents the contribution coming from the N elements in the double-sum, in which 
the action density in x and is evaluated in the field of the same instanton. One finds 

n , x _ M [dnidi(Sh) sjn^x) gg^o) 

where 

AT 

di(fii) = e- Gl ^) / ]Jdn 3 e-^ G ^ +G ^ ni \ (21) 
^ i>i 

The one-instanton distribution d(Q) cannot in general be computed explicitly from first 
orinciples, because of the breakdown of the perturbative expansion at large instanton sizes 



23j |. However, some of its general properties can be derived directly from symmetry consid- 
erations. In fact, translational invariance imply that d(Q) does not depend on the instanton 
position z. Moreover, gauge invariance of the action implies that does not depend on the 
color orientation. On the other hand, scale invariance is broken by the trace anomaly, hence 
d(Q) represents the instanton size distribution. 

In conclusion, the Hi contribution to the point-to-point correlation function can be ex- 
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pressed as 

N j d A z j dp di{p) s{p, x-z) s{p, -z) 
1[X) V fdpd(p) 

= (192) 2 r 7 U u x [ d 4 z [ dpd^p) P - t j, (22) 

1 > Jdpd(p)J J H lyH \( x ^ z )2 + p 2f ( z 2 +p2) 4' ^ > 

where n is the instanton density, which we shall not need to specify. 

Let us now consider the second contribution H 2 (x), which comes from all terms in the 
double sum in ( fl5l) . in which the scalar operator in x and in is evaluated in the field of 
different instantons. 

nix*., v AT(AT ^ /dgijgg d 2 (Q u Q 2 ) s(x;Q 1 ) g (0;f) 2 ) 
" 2 ( ' rl = 7V(7V L) /dfiidfiadCn^fla), 



N(N — 1 



/ d 4 2;idpidL r id 4 z 2 rfp 2 ^2 Pi, U u z 2 , p 2 , ?7 2 ) t( x -z^\ 



where 



/ d Z\dpidUid z 2 dp 2 dU 2 d 2 (zi, pi, Ui, z 2 , p 2 , U 2 ) 

(23) 



N 



d 2 (VL u Vt 2 ) = e - G i( n i)-Gi(n 2 )-G 2 (ni,n 2 ) /" d^e - ^> 2 G ^ n 0+ G 2(fh A)+ G 2(n 2 A) ^4) 

is the two-body instanton density. 

Some of the general features of the two-body density G? 2 (fii,fi 2 ) can be inferred from 
symmetry arguments. First of all, we note that the dependence on the relative color orien- 
tation of the instantons is in general non-trivial. However, such dependence does not affect 
the calculation of the action density correlator, because the two terms s(x; fli) and s(0; fl 2 ) 
are separately gauge invariant. Consequently the integral over the relative color orientation 
leads to an overall constant factor, which is canceled by the denominator. In addition, trans- 
lational invariance implies that the density does not depend on Z\ and z 2 separately, but 
only on the relative distance, z 2 — z\. Finally, we observe that infinitely separated instantons 
must be uncorrelated, so that the density correlation must approach a constant. Based on 
such considerations Il 2 can be written as 

/ ^ z x dp x A A z 2 dp 2 d' 2 (z! - z 2 , p u p 2 ) (( f 2 /, 1 2)2 ,^1% 

II 2 (z) = N(N- l) J - (( ° +P ; } ( 2+P2) , (25) 

J d ^idpid z 2 dp 2 d'(z 1 - z 2 , p 1 , p 2 ) 



where 



d' 2 (\z 1 -z 2 \, Pl , P2 )= I ' dUidu 2 d(nx,n 2 ). (26) 



n 



It is immediate to verify by inspection that the integrand is largest in the region in which 
\z\ — 22 1 ~ \x\. Hence, in the large \x\ limit, in which the two-body density becomes constant, 
the contribution il 2 (x) converges to 

n 2W OiZzll x (32)r) . (27) 

Hence, in the thermo dynamical limit N, V — > oo, we obtain a prediction for the square of 
the gluon condensate. 

In order to complete our calculation of the scalar correlator, we are left with specifying 
the one- and two-body density functions d\(p) and d' 2 {\z\ — z^\, pi, p-z). The instanton density 
d\{p) has been evaluated on the lattice, using different methods (see e.g. |30| and references 
therein). On the other hand, no lattice calculation for d' 2 {\z\ — z 2 |, pi, P2) has been performed 
so far. A possibility is to compute such a quantity in a model-dependent manner, for example 



using the choices for G\ and G2 of the Interacting Instanton Liquid Model 28| . which is 
briefly described in the appendix. In this work, however, we choose to avoid introducing 
additional model dependence and we try a different strategy, based on the Single Instanton 



Approximation (SIA) [26|]. The main idea is to evaluate the momentum-dependent effective 
mass in a regime of momenta p and Euclidean time r, in which the scalar correlator is 
dominated by Hi(x), hence it is sensitive mainly to the lattice-calculable function d\(p) and, 
not to the unknown two-body density d' 2 {\z\ — Z2I, Pi, P2) ■ 

IV. SCALAR GLUEBALL IN THE SINGLE INSTANTON APPROXIMATION 



The SIA was first developed in 12|, |26| and applied in [14I flo| to the calculation of electro- 
magnetic form factors. The starting point of this approach is to observe that, for sufficiently 
short-distances x, the scalar two-point correlation function U(x) is saturated by the single- 
instanton term n^x) — which depends on the lattice calculable distribution d\{p) — , and 
becomes insensitive on the multi-instanton term ^(x) — which depends on the unknown 
two-body density d 2 (Q,Q) — . 

The evaluation of the scalar glueball mass requires computing the momentum-projected 
correlation function ([3]) at large values of the Euclidean time r, where the contribution of the 
lowest bound-state can be isolated form from that of the higher scalar excitations. Clearly, 
as the Euclidean time comparable with the typical distance between two pseudo-particles 
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~ f?," 1 / 3 , we expect the point-to-point correlator H(x) — and therefore the momentum- 
projected correlator Gs(r, p) — to start receiving significant contributions from the two-body 
term n 2 (2;). 

The main idea of the present SIA calculation resides in observing that the range of 
Euclidean times in which the G$(p, t) is insensitive on the multi-instanton contribution 
Il2(x) can be increased, by projecting on states with large momenta p. 

With such a choice, one suppresses the dependence on the multi instanton density d' 2 (\zi — 
Z2\, Pi, p-i) which is responsible for the square of the gluon condensate. In addition, the 
momentum projection integral in ([3]) becomes sensitive only to the value of the point-to- 
point correlator II(x, r) in the spatial points at a distance < l/|p| from the Euclidean time 
axis. Hence, the effects of correlations with pseudo-particles which are at a distance larger 
than ~ l/|p| from such an axis is suppressed. 

The question we want to address in this section is whether there exist a window in p and 
r, in which the single-instanton term Hi(x) still dominates over the multi-instanton term 
n 2 (x), and the lowest scalar state dominates over its higher excitations. If this condition is 
verified, one should observe a range of Euclidean times in which the momentum-dependent 
effective mass M e //(p, r) defined in Eq. ([6]) is independent on both p and r. As it was 
argued in sect, flj], this would represent a clean signature of the existence of a scalar glueball 
bound-state in the instanton vacuum. 

Based on such considerations, we now compute the momentum projection correlator 
Gs{p,t) keeping only the one-instanton term n 1 (x). A number of independent lattice 
simulations have shown that the instanton distribution d\(p) has a peak in the range 
p ~ 0.35 — 0.45 fm, in good agreement with the estimate p ~ 1/3 fm, proposed long 
ago by Shuryak, based on phenomenological arguments Q. However, the details of the size 
distributions are quite sensitive to the prescription used to extract the instanton content 
from a generic vacuum gauge configuration. Hence, in this work, we choose to present the 
results obtained adopting the simplest form, 

d 1 {p) = n5{p-p), (28) 

where n ~ 3 fm -4 is the instanton density, which was extracted from Monte Carlo simulations 
in the Interacting Instanton Liquid Model 111 ] . We have checked that the results obtained 



using a smeared size distribution which interpolates one obtained in lattice calculations agree 
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within 10% with those obtained using the simple parameterization (128|) . Such differences 
are arguably smaller than the accuracy of the present calculation. 

After integrating over the instanton-size, the single-instanton contribution to the scalar 



point-to-point correlator reads 



21 



3l| 



1 1 



m st -(x) = (192) 2 np 8 / d 4 z j j. (29) 

3 1 ' 1 ' H J ((x-^+p 2 ) 4 (z 2 + p 2 ) 4 1 ' 

In order to construct the momentum-dependent effective mass ([6]), we take the spatial 

Fourier-transform of the calculated point-to-point correlator. The free contribution reads: 

1 fi n 2 

Gf e (r, |p|) = — f- e^l p '(3 + 3r|p| + r 2 p 2 ), (30) 

where the QCD coupling constant g has been determined at the chiral symmetry breaking 
scale, where the strong fine-structure constant is a s ~ 0.3. The instanton contribution in 
the SIA is 



Gf A (r, |p|) = 64vr 4 p 8 / dz,— — • 3 + 3 A 2 + p 2 |p| + V \z\ + p 2 ) 

J ((z 2 +p 2 )((r-Z 4 )2 + p2))V2 ^ V *->) 

■ (3 + 3^4 - r) 2 + p 2 |p| + p 2 ((^ 4 - r) 2 + p 2 )) . (31) 



From the momentum-projected correlation functions fl30l) and fl3Tl) it is immediate to 
construct the momentum-dependent effective mass: 



Meff(r,p) 



g(Gf-(r,p) + Gf A (T,p)) 
\\{ G| ee (r,p) + Gl IA (r,P) 



p 2 . (32) 



The results for M e //(r, p), calculated at different instanton sizes and different momenta 
are reported in Fig. [2] and Fig. [31 Fig. [2] presents the momentum-dependent effective mass, 
in a range of momenta, for few fixed instanton sizes. On the other hand, Fig. [2] shows the 
effective mass, in a range of instanton sizes, for few fixed momenta. These plots clearly show 
that there exists a range of Euclidean times in which the momentum-dependent effective 
mass becomes independent on both momentum and Euclidean time. As was argued in 
sect, mi this represents a signature of the existence of a scalar glueball bound-state in the 
instanton vacuum. 

On the other hand, at large Euclidean times and for very small momenta, we observe that 
the calculated effective mass starts depending on r and p again. This is well understood, 
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FIG. 2: Effective-mass plots evaluated at fixed instanton-size for different momenta. 
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FIG. 3: Effective-mass plots evaluated at fixed momentum, for different instanton-size. 

since in such regimes the contribution of a single instanton to the correlation function is no 
longer dominant and therefore the SIA breaks down. However, by stationary-phase analysis, 
we expect that if the value of the momentum p is increased, the integral ([3]) should become 
sensitive to shorter-sized point-to-point correlation functions. As a consequence, the range 
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FIG. 4: Scalar glueball mass calculated in the SIA, as a function of the instanton size p. 

of validity of the SIA should extend for longer Euclidean times. 

It should be stressed that, as the Euclidean time grows, we expect instanton-instanton 
correlations encoded in the contribution II 2 (x) to the two-point correlation function to be- 
come important. Since this part is neglected, the plateau in the analytic effective mass plots 
shown above eventually breaks down, at large Euclidean times. 

The complete list of masses extracted for different instanton sizes and different momenta 
are reported in Table I. We stress the fact that, for any given fixed instanton size, calculations 
performed at different momenta p lead to the same glueball mass. This represents an 
important consistency check of our calculation. 

In order to further assess the accuracy of our calculation, we have analyzed the dependence 
of the glueball mass on the the instanton size. Since the average instanton size is the only 
independent dimensional parameter in our model, the mass of the resulting scalar glueball 
should scale as 1/p. The plot in Fig.(j4j) shows that all our results can be fitted with very 
high accuracy by the function 

2 16 

M 0++ (p) = — - . (33) 

Our final prediction for the scalar glueball mass is obtained by choosing the instan- 
ton size in the range p ~ 0.30 — 0.33 fm, i.e. the value suggested by phenomenology 
and by lattice simulations. We find mass of M^++ ~ 1.290 — 1.420 GeV. Within its un- 
certainty, this result is consistent with the lattice result reported by Meyer and Teper, 
Mg+*+ = 1475(30)^.(65),,^. MeV. Such an agreement is quite remarkable, given the sim- 
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|p| = 2.45 


IpI = 4.9 


|p| = 7.35 


IpI = 9.8 


|p| = 12.25 


IpI = 14.7 


|p| = 17.15 


av. (err.) 


p = 0.04 














10.528 


10.528 (-) 


p = 0.08 










5.390 


5.392 


5.395 


5.392 (0.001) 


p = 0.12 




3.430 


3.535 


3.550 


3.553 


3.555 


3.555 


3.530 (0.020) 


p = 0.16 




2.707 


2.695 


2.707 


2.707 


2.710 


2.705 


2.705 (0.002) 


p = 0.20 


2.166 


2.166 


2.168 


2.171 


2.173 


2.171 


2.173 


2.171 (0.001) 


p = 0.24 


1.801 


1.798 


1.803 


1.803 


1.803 


1.806 


1.801 


1.803 (0.001) 


p = 0.28 


1.534 


1.548 


1.556 


1.556 


1.556 


1.566 




1.553 (0.005) 


p = 0.32 


1.333 


1.348 


1.350 


1.348 








1.345(0.002) 


p = 0.36 


1.186 


1.188 


1.181 


1.174 








1.183(0.002) 


p = 0.40 


1.068 


1.073 


1.058 










1.066(0.002) 


p = 0.44 


0.953 


0.956 












0.954(0.001) 


p = 0.48 


0.875 


0.848 












0.862(0.010) 


p = 0.52 


0.786 


0.750 












0.769(0.012) 



TABLE I: Glueball mass in GeV extracted at different probing-momentum (in GeV) and at different 
instanton-size (in fm). In all calculations the diluteness has been kept fixed to n ~ 0.03. The mass 
is only given for those values of |p| and p for which a plateau could be identified unambiguously. 



plicity of our calculation. It suggests that the binding of the compact scalar glueball mass 
is not strongly affected by confinement. We also note that our calculation agrees with the 
instanton-improved OPE prediction of Forkel M^ + E = 1.25(0.20) GeV. 

V. THE GLUEBALL IN THE INSTANTON VACUUM 

In this section, we comment on the results of our calculation and discuss at the qualitative 
level the dynamical mechanism underlying the binding of scalar glueballs, in the instanton 
vacuum. 

Our starting point consists in observing that the statement of the existence of a bound- 
state in the spectrum can be formulated as a statement about the dynamics of quantum 
de-localization of a point-like excitation, in the Euclidean space. In order to see this, we recall 
that a necessary and sufficient condition for the binding of a scalar glueball in gluondynamics 



17 



FIG. 5: The (normalized) scalar correlators in co-ordinate representation IT'(r, |x|), calculated at 
different Euclidean times: r = 0.1 fm (left panel), r = 0.5 fm (center panel) and r = 1.5 fm (right 
panel). The dashed lines correspond to the propagators in the free gluon theory, while the solid 
lines represent the propagator of a bound-state of mass M = 1.5 GeV. 

is that the scalar point-to-point correlation function behaves as 



where the first term is the gluon condensate, while the second term is the propagator of 
a scalar particle of mass Mq++ and \$ is the coupling of the scalar glueball state to the 
overlapping operator. 

The corresponding two-point function in momentum-space representation, defined in 
Eq. ([3]) decays like a single exponential at large Euclidean times, for any fixed momentum p, 



We note that this equation implies that the correlator Gs{t, p) will be dominated by lower 
and lower momentum components, as the Euclidean time increases. Hence, by uncertainty 
principle, Eq.s ( 155]) and (I3"4"j) express the progressive de-localization of the center-of-mass of 
a free single-particle quantum state. 

This behavior differs strongly from the characteristics of the free gluon theory. A gluonic 
excitation with ++ quantum numbers created in the origin at time r = 0, will diffuse 
according to the power-law 



M ++ 



= K 1 [M 0++ Vr 2 +^) 



(34) 



47T 2 \/t 2 + x : 




(35) 




(36) 



In momentum-space representation, this corresponds to the law 




(37) 
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The spatial dependence of the propagator in coordinate-space representation of a scalar 
excitation in the free theory and in a theory with a scalar glueball bound-state are confronted 
in FigJSJ We note that, for small Euclidean times, the free theory correlator is more peaked 
than in the case of bound-state correlator. On the other hand, at large Euclidean times, 
the situation is reversed: the free correlator is more spread-out than the single-particle 
correlator. In other words, as expected, free gluons de-localize more rapidly than a system 
of interacting gluons, bound in a glueball. 

Let us now argue on the picture which emerges from our instanton model calculation 
for gluon-dynamics. The vacuum of QCD and of Yang-Mills theory is characterized by the 
occurring of tunneling transitions, between degenerate classical vacua. During each of such 
transitions, gauge field fluctuations with solitonic structure emerge: the instantons. Since 
such fields are quasi-classical, they represent collective states made by a large number of 
gluons. On average, instanton fluctuations are isotropically distributed over space-time and 
hence give raise to a finite vacuum expectation value of the gluon field strength, i.e. the 
gluon condensate. At vanishing (four-) momentum, no further structure of the vacuum 
gauge field fluctuations is evident. 

At finite momenta, however, the local structures of the instanton fluctuations begin to 
manifest themselves. For sufficiently short times, the short-distance correlations introduced 
by a single-instanton — Eq. ff29l) — dominate over the power-law correlations of the free- 
theory — Eq. (|36l) — . As shown in Eq. (I3li) . the occurrence of such a tunneling transition 
pumps-in additional high-momentum components into the momentum-projected correlator 
(1371) . Our results have shown that such components are sufficient to locally turn the free- 
gluon de-localization law — Eq. (1371) — into a single-particle de-localization law — Eq. (|35|) — . 

Such a dynamical correlation can be pictured as following: during a quantum vacuum 
rearrangement, gluons interact by exchanging momentum with the field of an instanton. 
Obviously, a single tunneling event is not sufficient to correlate the field of a glueball for 
arbitrarily large Euclidean distances. In order to generate a state which diffuses like a 
single-particle state, for arbitrary Euclidean times and arbitrarily low-momenta, one needs 
to re-sum the contribution of an infinite sequence of tunneling transitions, i.e. to account 
for an infinite number of scattering processes of the gluons with the instanton fields. 
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VI. CONCLUSIONS AND OUTLOOK 



In this work, we have studied the short-distance instanton-induced interactions in the 
scalar glueball channel. In particular, we have addressed the question whether such inter- 
actions are sufficient to generate a glueball bound-state with ++ quantum numbers and 
realistic mass. 

We performed a SIA calculation in the quenched approximation, which in this channel 
corresponds to gluon-dynamics. We observed evidence for a bound-state and we have esti- 
mated its mass to be M^++ ~ 1.290 — 1.420 GeV, i.e. consistent with the value calculated 
by Meyer and Teper in lattice gauge theory. We stress the fact that our results do not rely 
on any a priori assumption on the structure of the spectral function. The emergence of a 
plateau in the effective mass plot in this model represents a strong indication that the origin 
of the scalar glueball is mostly provided by short-sized non-perturbative gluon fluctuations, 
rather than by confining forces. 
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APPENDIX A: INTERACTING INSTANTON LIQUID MODEL 



In the following, we shall illustrate a particular form of the one- and two-body correla- 
tions Gi(Q\) and Gr 2 (f2i, uj 2 ), which appear in Eq. (jTTI) . This choice defines the so-called 



Interacting Instanton Liquid Model 



28] 



The one-body term Gi(Q), is derived from the quantum fluctuations around the single 
instanton configuration is calculated in Gaussian approximation, through 't Hooft's 
semi-classical instanton amplitude d(pi) 23] 

b'\ b' 



G\ oc — In 



P~V V ^i(pr Vc exp -P 2 (p) +[2N C + 



2b J 2b(3 1 (p) 



(Al) 



The functions j3i and (5 2 are defined as 



Pi(j>) = -bHpA PV ) p 2 (p) = + ^ln (lfr{p) 
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where Apy = 1.03 is the QCD scale in the Pauli-Villars scheme. The coefficients 
6, V are given by b = ^N c - |JV> and V = f N c 2 - f JV c JVy + J. 

A classical two-body interaction between instanton and anti-instanton G^' is derived 
from the Yang-Mills action evaluated in the stream-line ansatz (jA3j) . 

G£{n u Qj) = S YM [Af eam in h - 2S , (A2) 

where So = ^ is the action of an isolated instanton. The total stream-line gauge field 
is constructed by summing-up the field of instanton-antiinstanton pairs. The analytic 
gauge potential for a single instanton-antiinstanton pair is given by 



A "^ 2 ^^x +2R ^4 ^Tm < A3) 



where A is the streamline conformal parameter defined by 



R 2 + p 2 i + p 2 a , ( (R 2 + p 2 i + p 2 a) 2 a 1/2 

ZpipA V ±p\p\ 

R = \zi — za 



a = AL n J „ A + r \ h Ft - 1 (A4) 



and 



R ab = -Tr[U^r a UT b ] (A5) 
represents the relative color orientation of the two instantons. 

The effect of quantum two-body correlations is accounted for by means of a phe- 
nomenological short-distance hardcore term [28| : 

Gf e "-(a,^) = ^M 2 , (A6) 

where u is a color orientational factor. Such a term removes large-sized instantons 
from the vacuum and provides a cut-off to the momentum that can be exchanged 
through the instanton field. Hence, it restricts the region of applicability of the model 
to the non-perturbative sector, characterized by momenta of the order p < 1/p, where 
p is the average instanton size. The dimensionless coefficient A in (IA6I) controls the 
strength of the repulsion and is the only phenomenological parameter of the model. 



Most calculations adopt the value A ~ 128, first suggested by Schafer and Shuryak [28] 



and used also in 
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